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Abstract 

We show that topological frequency band structures emerge in two-dimensional electromagnetic 
lattices of metamaterial components without the application of an external magnetic field. The 
topological nature of the band structure manifests itself by the occurrence of exceptional points in 
the band structure or by the emergence of one-way guided modes. Based on an EM network with 
nearly flat frequency bands of nontrivial topology, we propose a coupled-cavity lattice made of 
superconducting transmission lines and cavity QED components which is described by the Janes- 
Cummings-Hubbard model and can serve as simulator of the fractional quantum Hall effect. 

PACS numbers: 42.70.Qs, 73.20.-r, 73.43.-f 
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I. INTRODUCTION 



The topological description of the quantum states of matter sets in a new paradigm 
in the description and classification of atomic solids. Namely, atomic solids whose energy 
band structure possess nontrivial topological properties constitute a new class of materials 
whose salient properties are robust to phase transitions which modify the symmetry order 
of the atomic solid. Prominent examples of such topological atomic solids are the integer/ 
fractional quantum Hall (I/FQHE) systems and the topological insulators (TIs). Well-known 
examples of topological properties are the existence of chiral edge states in QHE systems 
and the presence of gapless surface states in TIs which are both immune to order-disorder 
phase transitions. 

The advent of artificial electromagnetic (EM) structures such as photonic crystals and 
metamaterials has established over the years a continuous conveyance of ideas and methods 
from atomic solids to their EM counterparts. Quite naturally, the concept of topological 
order has been adapted to photonic crystals starting with the QHE: a two-dimensional 
(2D) lattice of gyromagnetic/ gyroelectric cylinders is a system with broken time- reversal 
symmetry^ with frequency bands characterized by a nonzero Chern number, allowing for the 
emergence of unidirectional (one-way) edge states^ in analogy with the chiral edges states 
in QHE systems such as a 2D electron gas or graphene nanoribbons under magnetic field. 
Anomalous QHE can also be simulated with artificial chiral metamaterials of gyromagnetic 
components.- In TIs^ and quantum spin Hall systems 5 in 2D, the presence of magnetic field is 
not prerequisite for the appearance of topological electron states. In analogy with atomic TIs, 
in certain 3D photonic crystals and metamaterials with proper design, topological frequency 
bands appear without comprising gyromagnetic/ gyroelectric materials which require the 
application of external magnetic field in order to break time-reversal symmetry.- - - 

In this Letter, we propose a class of 2D EM networks possessing topological frequency 
bands without the application of an external magnetic field. Namely, we show that topo- 
logical bands emerge in 2D lattices of EM resonators connected with left- and right-handed 
metamaterial elements such as transmission lines or waveguides loaded with a negative 
refractive- index medium. The topological nature of the corresponding frequency bands is 
manifested by the emergence of an exceptional point for transverse electric (TE) and by the 
generation of one-way modes for transverse magnetic (TM) waves. In the latter case, the 
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system can be viewed as a simulator of the FQHE for polaritons. 

II. LATTICE OF COUPLED DIPOLES. 

The EM crystals under study here are amenable to a photonic tight-binding description 
within the framework of the coupled-dipole method.^ The latter is an exact means of solving 
Maxwell's equations in the presence of nonmagnetic scatterers. We consider a lattice of 
cavities within a lossless metallic host. The i-th cavity is represented by a dipole of moment 
Pj = (Pi- X , Pi-y, Pi-z) which stems from an incident electric field E mc and the field which is 
scattered by all the other cavities of the lattice. This way the dipole moments of all the 
cavities are coupled to each other and to the external field leading to the coupled-dipole 
equation 



Gn'(u) is the electric part of the free-space Green's tensor and oti(uj) is the polarizability 
the i-th cavity. Eq. ([T]) is a 3N x 3N linear system of equations where N is the number of 
cavities of the system. 

For a particle/cavity of electric permittivity e embedded within a material host of 
permittivity e^, the polarizability a is provided by the Clausius-Mossotti formula a = 
(3V/47r)(e — e h )/(e + 2e h ), where V is the volume of the particle/ cavity. For a lossless 
Drude-type (metallic) host i.e., €h(co) = 1 — Up/u 2 (where u p is the bulk plasma frequency), 
the polarizability a exhibits a pole at lu = uo v ^/2j[t + 2) (surface plasmon resonance). By 
making a Laurent expansion of a around uq and keeping the leading term,- we may write 
a = F/(u - u Q ) = F/tt where F = (27V '/8n)tu e/ '(2e + 4). For sufficiently high value of the 
permittivity of the dielectric cavity the electric field of the surface plasmon is much localized 
at the surface of the cavity. As a result, in a periodic lattice of cavities, the interaction of 
neighboring surface plasmons is very weak leading to much narrow frequency bands. By 
treating such a lattice in a tight binding (TB) manner,- we may assume that the Green's 
tensor Gu'(oj) does not vary much with frequency and therefore, Gtf(u) ~ Ga'(uo). In this 
case, Eq. ([!]) becomes an eigenvalue problem 




(1) 




(2) 
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where F has been absorbed within the definition of Gu>(ujq) an d we have set E mc = in 
Eq. (JT]) as we are seeking the eigenmodes of the system of cavities. In the following, we will 
be dealing with 2D lattices of cavities. We can, therefore, treat separately the case where 
the electric field lies within the plane of cavities (TE modes) from the case where the electric 
field is perpendicular to the plane (TM modes). 



A. TE modes 



In this case, Pj = (Pi- X , Pi ;y ) and the Green's tensor Gu>(uo) is given by 
Gtf/(w ) = Fql C(<7oN'|)l2 + J{qo\ru>\) 



x ii'Vii' y ii> 



(3) 



with r&/ = rj — iv, q = Wovft(wo)/c and I2 is the 2x2 unit matrix. Since we focus our 
attention around the surface plasmon frequency ojq, we operate in the subwavelength regime 
where (faKi'l -C 1. In this regime, the functions C(q \ra>\), J(qo\ra>\) are written as 

Qo FC (Q\\\ r u'\) - -9o^(9||kii'|) - q% F exp(iqo\rn>\)/(qo\ru>\) 

= ta> exp(i<fi u ,) (4) 

where tu> and are real numbers. In what follows, the cavities are connected via coupling 
elements, i.e., waveguides or transmission lines, in which case the phase factors 0y are 
not necessarily related with the wavevector of the host medium and can therefore be 
considered as independent parameters. 

For a 2D lattice of cavities, we assume the Bloch ansatz for the polarization field, i.e., 



n/3 



exp(ik • RrjPo/? (5) 



The cavity index i becomes composite, i = n/3, where n enumerates the unit cell and j3 the 
positions of inequivalent cavities in the unit cell. Also, R„ denotes the lattice vectors and 
k = (k x , k y ) is the Bloch wavevector. By substituting Eq. (jHJ) into Eq. (J2J) we finally obtain 

^G j g /9 /(wo,k)P ^ = nPo^ (6) 

0' 



4 



where 

G/3/j'(u;o,k) = ^exp[ik ■ (R n - Rn')]G n p in rpr(u ). (7) 

Solution of Eq. ([6]) provides the TE frequency band structure of a periodic system of cavities. 

In order to seek for topological Bloch modes in a 2D lattice, we need at least two distinct 
frequency bands. Since, TE modes correspond to two degrees of freedom for the polarization 
field, i.e., (P x ,P y ) we may consider a 2D lattice with one cavity per unit cell. Namely, we 
consider the square lattice of Fig. [Th^ where we consider nearest-neighbor (NN) and next- 
nearest-neighbor (NNN) hoppings of the EM field among the cavities. The NN hopping 
carries a nonzero phase texp(±0) whose signs are denoted by the arrows in Fig. [Th. The 
NNN hopping is denoted by if. For the lattice of Fig. [T^,, the Green's tensor G of Eq. (jj[J) 
becomes 

G = t[cos(f>cos(k x a/2) cos{k y a/2) — ism.(j)sm.{k x a/2) sin(A; 2/ a/2)] 
\ y ) 

cos(k x a) I 

Fig. [Tb shows the frequency band structure derived from Eq. (jSJ) for t — t' — 1, <p = 
We observe that at some point along the XM symmetry line the frequency bands coalesce 
into a single band. The point beyond which the bands coalesce is an exceptional point 
and has been observed in PT-symmetric lattices.— In general, exceptional points emerging 
in parameter space are associated with topological charge and geometric (Berry) phase.— 
The topological properties of an exceptional point have been revealed by encircling it in 
parameter spaced as it was demonstrated in a microwave cavity experiment.— Although a 
proper theory of the topological properties of the exceptional points in lattices is still lacking, 
based on previous work 11 we can indirectly assign topological charge and geometric phase 
to the exceptional point appearing in the frequency band structure of Fig. [Tb. 

B. TM modes 

Next, we assume that the polarization at each dipole is oriented in the z-axis. In this 
case, the eigenvalue problem of Eq. (JSJ) becomes scalar and 

G& = FqlC{q\\\r iV \) ~ go^ ex P(^ohd)/(<?oh;'|) = tw exp(i<f> & '). (9) 
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FIG. 1: (Color online) (a) Square lattice of EM resonators connected with metamaterial-based 
coupling elements. The arrows denote nearest-neighbor hoppings whilst the solid lines next-nearest- 
neighbor hoppings. The direction of the arrow shows whether wave propagation in the coupling 
element is left- or right-handed, (b) Frequency band structure corresponding to lattice of the left 
panel for t = t' = 1, <fi = 7r/3. 

The same applies to Eqs. ©-© and the TM problem becomes equivalent to the electronic 
case. Since the minimal model to have topological frequency bands is a two-band model, 



we adopt the checkerboard lattice of Ref. [14| (see Fig. (2^,). Namely, apart from considering 



NN and NNN hoppings as in Fig. we also consider next-next-nearest-neighbor (NNNN) 
hoppings (denoted by the arcs in Fig. [2^) with strength t" . The NN hoppings are, again, 
complex, texp(±0) where the sign is denoted by the arrows in Fig. The NNN hopping 
strength is t\ (t' 2 ) if two sites are connected by a solid (dashed) line. We note that in EM 
lattices such as those considered here, a negative phase — can be easily achieved when the 
cavities are connected e.g., by ID left-handed transmission lines (LHTL), i.e. transmission 
lines supporting backward-propagating waves where the phase velocity is opposite to the 
group velocity.— 1 ^ Alternatively, the cavities may be connected by waveguides loaded with 
a left-handed (LH) metamaterial. Obviously, a positive phase +<ft can be achieved by similar 
means [right-handed transmission lines (RHTLs)]. 

For the checkerboard lattice of Fig. [2^, the Green's tensor Gbb 1 becomes 



G = (10) 
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FIG. 2: (Color online) (a) Checkerboard lattice of EM resonators connected with metamaterial- 
based coupling elements. The arrows denote nearest-neighbor hoppings whilst the solid and broken 
lines next-nearest-neighbor hoppings. The direction of the arrow shows whether wave propagation 
in the coupling element is left- or right-handed. Two of the next-next-nearest neighbor hoppings 
are shown as dotted arcs, (b) The EM resonator of the lattice is a superconducting circuit QED 
system consisting of an LC resonator coupled to Cooper pair box (CPB). Typical transmission line 
for (c) right- ((d) left-) handed coupling elements along with the corresponding dispersion relation. 

where 



Gii = 2ti cos{k x a) + 2t' 2 cos(k y a) + 4t" cos(k x a) cos(k y a) 
G\2 = G* 2l = 4tcos(ficos(k x a/2) cos(k y a/2) — 4it sin <psxn(k x a/ 2) 
G 2 2 = 2t' 2 cos(k x a) + 2t[ cos(k y a) + At" cos(k x a) cos(k y a) 



(11) 
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Atomic lattices with exotic hoppings such as those considered here have been used for 
simulating fractional Quantum Hall effect (FQHE) states at zero magnetic field as nearly 
flat topological bands can emerge which simulate the Landau levels associated with a uniform 
magnetic field.— _ — By taking as TB parameters,— t — 1, = 7r/4, t' x — —t' 2 = 1/(2 + \/2), 
t" = 1/(2 + 2y/2)], a nearly flat band emerges as it is evident from the frequency band 
structure of Fig. [3k. Based on the equivalence of the Green's tensor of Eqs. (flQ|) and (|TT|) 
with the Hamiltonian of the electronic problem,~ each of the two bands of Fig. [3^ carries a 
Chern number ±1. The topological nature of the frequency bands of the lattice of Fig. [2h, is 
also manifested by the emergence of one-way bands (the photonic counterpart of the electron 
chiral edge states^) in the frequency band structure for a slab geometry of Fig. [3b. 
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FIG. 3: (a) Frequency band structure for the infinite checkerboard lattice of Fig. [2] and (b) the 
frequency band structure for a finite slab consisting 20 unit planes [parameters: t = 1, <p = tt/4, 
t[ = -t' 2 = 1/(2 + V2), t" = 1/(2 + 2v/2)] 

The occurrence of topological properties such as the exceptional point in Fig. [Tb and the 
one-way modes in Fig. [3] are a result the synthetic gauge field which is generated by the 
geometry of the metamaterial-based coupling elements (formation of closed flux loops of the 
phase of the EM field). 
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III. SIMULATION OF THE FQHE 



Having established a nearly flat topological frequency band of the EM field for the lattice 
of Fig. [2^, we are able to design a system for creating an EM analog of the FQHE. The 
most natural choice would be to consider a coupled cavity array (CCW) wherein polaritons 
propagate through a hopping mechanism (as in our case) and interact strongly with the 
reservoir of modes when they reside within the cavity.— >2i FQHE with magnetic field can 
also be simulated by atoms confined in a 2D CCW. 22 As stated above, the EM lattice of 
Fig. [2^ can be realized in the laboratory as a network of transmission lines (TLs) where 
hoppings with positive (negative) phase can be realized with RHTLs (LHTLs) as shown in 
Fig. [2t (Fig. |2H). This means that the topological bands lie in the GHz regime. Therefore, 
in order to simulate the FQHE for microwave photons we need to implement a cavity QED 
scheme in this regime. This can be achieved by considering a superconducting-circuit cavity 
QED system 2 ^ 2 ^ consisting of a Cooper-pair box (CPB) coupled to a TL resonator (see the 
equivalent circuit of Fig. |2b). The CPB operates as an artificial atom (two-level system )2£i2& 
and couples to the microwave photons of a superconducting TL resonator which plays the 
role of an on-chip cavity reservoir. The microwave response of superconducting circuit cavity 
QED system is described by the Jaynes-Cummings Hamiltoniar>22>2i 

H JC = hwriafa, + 1/2) - ~(£ e ^ + Eja\) (12) 
+ hg{a\or + ai af) 

where u r = 1/yLC is the frequency of the superconducting resonator, af (a^) creates (an- 
nihilates) a microwave photon in the TL resonator (cavity), af (a.;) creates (annihilates) an 
excitation in the CPB, g is the coupling parameter between the CPB and the TL resonator, 
E e i is the electrostatic energy and Ej = Ej )Tnax cos(7r$;,) the Josephson energy of the CPB. 
$5 = $/<l>o is a flux bias applied by a coil to the CPB and controls the Josephson energy 
Ej. 

A superconducting circuit cavity QED system where microwave photons propagate in the 
lattice of Fig. [2^ is described by a Jaynes-Cummings-Hubbard Hamiltonian of the form 2 ^ 

h jch = Y j Hf + H TB (13) 

i 

where H TB is the tight-binding form of the Hamiltonian of the microwave photons propa- 
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gating within the checkerboard lattice of Fig. [2^, i.e., 



H 



TB 



t y] exp (ifaj) {aj~ a j + H.c.) - 2J t 'ij{ a t a i + H.c.) 

(hi) {(hi)) 

t" ( a t a j + H.c.) 



(14) 



«<M»> 



which is the direct-space representation of the Green's tensor of Eq. ([TP]) . An important 
ingredient which gives rise to the FQHE is the presence of repulsive interactions among the 
microwave photons and is inherently present in Eq. (|13p as photon blockade.— The latter 
phenomenon has been recently observed experimentally in the GHz regime for superconduct- 
ing circuit cavity QED systems such as the one considered here (CPB + TL resonator).—^ 
The different FQHE phases can be calculated by direct-diagonalization of the Hamiltonian 
of Eq. (fl3l) . We note that the proposed quantum simulator for the FQHE differs fundamen- 
tally with previous proposals^^i since it essentially constitutes a passive design requiring 
no externally applied electric or magnetic fields. 

Some typical values of the Hamiltonian of Eq. (TX5T) are: 23 u r = 38 GHz, Ej tTnax = 8GHz, 
Eq = 5.2GHz, g 0.314 GHz. The latter parameter, g, is much larger than the loss rate of 
the TL resonator (~ 0.005 GHz) and the decoherence rate of the CPB (~ 0.004 GHz). The 
frequency u r of the TL resonator should fall within the operating bandwidth of the LH- and 
RHTLs. In Fig. [2b and Fig. [2H we have considered ideal TL which have infinite bandwidth. 
However, actual LH- and RHTLs have very large bandwidth which is a distinctive feature of 
nonresonant metamaterials compared to the resonant ones. 16 Lastly, if the superconducting 
QED chip has a thickness of 1mm, the coupling TLs (RH or LH) have 10mm length and 
the TL resonator covers an area of 30mm 2 , for the given resonator frequency (38GHz), the 
hopping strength t is about 0.5GHz which is also significantly larger than both the TL loss 
and CPB decoherence rates. 

In order to probe experimentally the FQHE with the proposed structure one needs to 
create the phase diagram of the spectrum gap between the FQHE ground-state manifold 
and the lowest excited states, as a function of the coupling parameters g for NN and NNN 
hopping when the latter lie in the photon blockade regime. Generally speaking, in the FQHE 
state the spectral gap assumes much larger values than in superfluid and solid phases.— 8 The 
frequencies of the ground and excited states (and thus their corresponding gaps) can be 
measured by microwave transmission experiments. 
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IV. CONCLUSION 



In conclusion, we have shown that topological frequency bands emerge in 2D electromag- 
netic lattices of metamaterial components in the absence of an applied magnetic field. The 
topological nature of the corresponding band structures gives rise to significant phenomena 
such as one-way waveguiding and coalescence of EM modes. The above lattices can be the 
basis for realizing a simulator for the FQHE based on superconducting transmission lines 
and circuit cavity QED systems. 



Electronic address: vyannop@upatras.gr 
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